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ABSTRACT 

We have carried out simulations to predict the performance of a new space- 
based telescopic survey operating at thermal infrared wavelengths that seeks 
to discover and characterize a large fraction of the potentially hazardous near- 
Earth asteroid (NEA) population. Two potential architectures for the survey 
were considered: one located at the Earth-Sun LI Lagrange point, and one in a 
Venus-trailing orbit. A sample cadence was formulated and tested, allowing for 
the self-follow-up necessary for objects discovered in the daytime sky on Earth. 
Synthetic populations of NEAs with sizes as small as 140 m in effective spher¬ 
ical diameter were simulated using recent determinations of their physical and 
orbital properties. Estimates of the instrumental sensitivity, integration times, 
and slew speeds were included for both architectures assuming the properties of 
newly developed large-format 10 /rm HgGdTe detector arrays capable of operat¬ 
ing at ~35 K. Our simulation included the creation of a preliminary version of a 
moving object processing pipeline suitable for operating on the trial cadence. We 
tested this pipeline on a simulated sky populated with astrophysical sources such 
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as stars and galaxies extrapolated from Spitzer Space Telescope and Wide-field 
Infrared Explorer data, the catalog of known minor planets (including Main Belt 
asteroids, comets, Jovian Trojans, planets, etc.), and the synthetic NEA model. 
Trial orbits were computed for simulated position-time pairs extracted from the 
synthetic surveys to verify that the tested cadence would result in orbits suitable 
for recovering objects at a later time. Our results indicate that the Earth-Sun 
LI and Venus-trailing surveys achieve similar levels of integral completeness for 
potentially hazardous asteroids larger than 140 m; placing the telescope in an in¬ 
terior orbit does not yield an improvement in discovery rates. This work serves as 
a necessary first step for the detailed planning of a next-generation NEA survey. 


1. Introduction 


From the kilometer-diameter Meteor Crater to the enormous ring-shaped feature be¬ 


neath the Gulf of Mexico resulting from the Cretaceous-Paleogene impactor( Alvarez et ah 


1980 Hildebrand et al.||1991 ), numerous lines of evidence record the history of previous near- 


Earth object (NEO) impacts. More recent events, such as comet Shoemaker-Levy 9’s impact 
with Jupiter (Levy et al. 1995), have spurred international interest in observational efforts 
to determine the collisional hazard to Earth posed by NEOs. To date, most efforts have 
focused on ground-based surveys carried out by 0.5 to 2 m class telescopes at visible wave¬ 
lengths (e.g. the Near Earth Asteroid Tracking Program, the Lincoln Near-Earth Asteroid 
Research Program, the Lowell Near-Earth Object Survey, Spacewatch, Pan-STARRS, and 
the Catalina Sky Survey; Helin et al.||1997 Stokes et al. 2000 Koehn fc Bowell||2000 ; Larson 


2007 McMillan 2007). These successful programs have now discovered 90% of all NEOs 
larger than 1 km in diameter (Mainzer et al. 2011b). A total of ~12,000 NEOs have been 
found to date {http://neo.jpl.nasa.gov/stats.html). Dedicated systems at the International 
Astronomical Union’s Minor Planet Center (MPC) archive and link new observations of 
NEOs in real time. Impact probabilities are updated in real time by the NASA Near Earth 
Object Program’s SENTRY system and the European Union’s NEODys system]^ ([Milani 


et al. 2005). 


The ~ 12,000 NEOs known so far are believed to represent only a small fraction of the 
total population. Statistical approaches to quantifying the hazard from the undiscovered 
population indicate that while the risk of a very large impact is low over the next several 
thousand years, dangerous NEO impacts remain a stochastic process on human timescales. 


^http://neo.jpl.nasa.gov/risk /; http://newton.dm.unipi.it/neodys/ 
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In 2004, the U. S. Congress mandated that NASA should hnd more than 90% of all NEOs 
larger than 140 m in diameter by the year 2020. The process of surveying for NEOs has 
proven extremely beneficial to planetary science. The survey data used to discover, catalog 
and characterize NEOs are responsible for finding the majority of comets. Main Belt aster¬ 
oids, Jovian Trojans, and Centaurs known today, a tremendous boon to those interested in 
understanding the origins and evolution of our solar system. 


Additionally, the NEOs that are likely to be the most hazardous tend to make numerous 
close approaches to the Earth, making some of them potentially attractive destinations for 
future human and robotic missions. NASA has been charged by the President of the United 
States with a new initiative to send human explorers to a near-Earth asteroid in the 2025 
timeframe (c.f. Abell et ah 2009). Out of the NEOs known to date, only a handful meet 


requirements for size and accessibility in the timeframe of interest with sufficiently short 
trip durations (http://neo.jpl.nasa.gov/nhats/). Little is known about these few objects 
in most cases other than their orbits, which are often highly uncertain, and their absolute 
visible magnitudes. NEO surveys are needed to find and characterize suitable targets for 
exploration, to identify and characterize potentially hazardous objects, and to carry out 
larger-scale scientific studies of the small body population of our solar system. All of these 
reasons motivate efforts to carry out a survey capable of not only discovering a large fraction 
of NEOs big enough to cause major damage but also to systematically characterize their 
physical and orbital properties. 


There is precedent for using space-based IR telescopes to discover and characterize 
asteroids in large numbers. NEOs emit a significant fraction of their bolometric luminosity 
at thermal infrared (IR) wavelengths, and the stellar and galactic backgrounds against which 
they are detected are typically much dimmer than in visible light. NEO thermal emission is 
only a weak function of an object’s geometric albedo, unlike its visible light brightness, so a 
telescope surveying at IR wavelengths is essentially equally sensitive to low and high albedo 
objects ( Mainzer et al.|2011b ). Moreover, surveying at IR wavelengths offers the opportunity 
to constrain an object’s diameter with much improved accuracy over visible wavelengths 
alone. Independent IR surveys are therefore useful when extrapolating the observed sample 
to determine properties of the population such as size frequency distributions. No conversion 
between H and diameter is needed, since diameter is determined directly, eliminating the 
uncertainty associated with the large range of possible asteroid albedos. If both visible and 


IR fluxes are obtained, it is possible to compute geometric albedo (Harris 1998; Mainzer 


et ah 2011cp(). Finally, IR wavelengths can be used to set constraints on the properties of 


regolith such as thermal inertia and conductivity (e.g. Ah-Lagoa et ah |2014; Groussin et ah 


2011 Delbo et ah 2007). 
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NEOWISE, the asteroid-hunting portion of NASA’s Wide-field Infrared Survey Explorer 


(WISE; Wright et ah 2010) mission, discovered and characterized >158,000 asteroids and 


comets throughout our solar system during its year-long prime mission (Mainzer et ah 2011a); 
>34,000 were new discoveries. Included in the sample were ~700 NEOs, of which 135 were 
new discoveries (Mainzer et al.i2011b, 2012a, Mainzer et ah|20i4 ). NEOWISE discovered the 
first known Earth Trojan, 2010 TK 7 ( Connors et ah||2011 ) as well as a so-called “horseshoe” 
NEO (2010 SOie; Christou fc Asherj2011[ ). By virtue of its wavelengths, all-sky coverage, and 
uniformity, the NEOWISE sample has been used to set constraints on the numbers, sizes, 
and orbits of asteroids, including NEOs and those that are classified as potentially hazardous 


(Mainzer et ah 2012bi Grav et ah 2011b, 2012). NEOWISE data have been used to study 


the size and albedo distributions of asteroids throughout the Main Belt (e.g. Masiero et ah 


2011 DeMeo fc Carry||2014), and to identify new collisional family members using albedo as 


an additional constraint (Masiero et ah 2013; Carruba et ah 2013 Walsh et ah 2013; Milan! 


et ah 2014). After a 32 month hibernation following the completion of its prime mission in 


February 2011, the spacecraft was brought back into operation, renamed NEOWISE, and has 
resumed discovering and characterizing NEOs using its remaining 3.4 and 4.6 /rm channels 


(Mainzer et ah 2014). 


Other space missions have been used to study asteroids and comets, with smaller num¬ 


bers of new discoveries. The Infrared Astronomical Satellite (IRAS; Neugebauer et ah 1984 


Beichman et al.||1988) surveyed the sky and observed ~2000 small bodies, mostly previously 


known Main Belt asteroids (Tedesco et al. 

2002 b 

). Additional observations came from the 

Mid-Course Space Experiment (MSX; Mill et al. 

1994 

Price et al.|2001 

Tedesco et al.|2002a). 

The Spitzer Space Telecope has observed hundreds of previously known NEOs (e.g. 

Trilling 


et al. 2010 Mueller et ah 2011; Thomas et ah 2011). The AKARI mission observed ~5000 


previously known asteroids (Usui et al. 2011, 2013 Hasegawa et al. 2013). 


Using a space-based telescope operating at thermal IR wavelengths to discover and 
characterize large numbers of NEOs is not a new concept 


Tedesco et al. 

... 

), 

Cellino 


et al. (2000), Cellino (2003), and Cellino et al. (2004) considered space-based thermal IR 


telescopes observing from Earth orbit, the Earth-Sun L2 Lagrange point, and a Venus-like 
orbit using HgCdTe detectors cooled by mechanical cryocoolers. However, the 2003 report of 


the NASA NEO Science Dehnition Team (Stokes et al. 2003) opted not to consider a space- 


based IR system on the grounds that sufficiently large, high-operability long-wavelength focal 
plane arrays did not exist at the time. At the time, the largest long-wavelength (~ lO/im) 
HgCdTe detectors designed for low-background astronomical background applications were 
in a 512x512 format, and their dark current and operability were inadequate for natural 
background-limited performance. Nevertheless, the initial results of that pilot program lead 
by the University of Rochester and Teledyne Imaging Systems were encouraging ([Bacon et al. 
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Motivated by the possibility of discovering NEOs in thermal IR wavelengths with WISE 
but aware of its limitations in terms of lifetime and held of view, our group revisited the space- 
based NEO survey designs. The major advance that resulted in WISE’s improved sensitivity 
and spatial resolution over IRAS was the increase in its focal plane formats to 1024x1024 
from 62 pixels, despite WISE’s smaller primary mirror (40 cm vs. 60 cm). Moreover, the 
focal planes of the Spitzer Space Telescope ( Werner et al.||2004 ) were predicted to equilibrate 
at ~30 K after the liquid helium was exhausted, a prediction that was borne out when 


the focal planes stabilized at 29 K after the cryogen ran out in 2008 (Storrie-Lombardi & 


Dodd 2012). This prediction suggested that if 10 pm HgCdTe arrays could be enlarged to 
a 1024x1024 format and could achieve high operability at ~35 K, it would be possible to 
cool the detectors and telescope purely passively. An orbit that maintains enough separation 
from the Earth to minimize its heat load while still maintaining a nearly constant distance, 
thus supporting high-bandwidth communications, can be found at either of the Earth-Sun 
LI or L2 Lagrange points. While WISE is constrained by its orbit and sunshade to viewing 
a narrow strip close to 90° solar elongation, a mission at LI with a taller sunshade can point 
much closer to the Sun. 

The resulting mission concept, called the Near-Earth Object Camera (NEOCam) was 
hrst proposed to NASA’s Discovery program in 2005 and again in 2010. In 2010, the project 
was awarded technology development funding to mature the long-wavelength IR HgCdTe 
arrays. The results of that effort have produced new 10 pm-cutoff 1024^ arrays with very 
high pixel operability at 35-40 K ( McMurtry et al.|2013 ). The mission design allows for long 
lifetime, high data transfer rates supporting downlink of full-frame images, and the ability 
to view large swaths of Earth’s orbit instantaneously. By enabling full-frame downlinks, 
standard astronomical data processing techniques for extracting sources and producing ac¬ 
curately calibrated astrometry and photometry can be employed. Many objects are detected 
at or near the detection threshold, so the ability to extract faint sources and distinguish 
them from noise and artifacts is essential. With NEOWISE, we have demonstrated the abil¬ 
ity to extract sources down to signal-to-noise (SNR) of 4.5 with extremely high reliability 


and completeness (Cutri et al. 2012 Mainzer et ah 2011a). Similar science data processing 


techniques can be used for NEOCam. 

Alternative orbits for carrying out an advanced NEO survey have been considered. For 


example, Tedesco et al. (2000) and Cellino et ah (2004) considered a Venus-like orbit on the 


grounds that the shorter synodic period with respect to Earth might afford more observations 
of NEOs with particularly Earth-like orbits. In theory, such an orbit might allow NEOs to 
be discovered more rapidly than a system operating at LI. However, a Venus-like orbit takes 
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a spacecraft ~30-170 times further from the Earth relative to LI. 

At present, Ka-band (the current state-of-the-art for spacecraft telecommunications 


systems, as demonstrated by the Kepler mission; Koch et ah 2004) can support ~150 Mbps 
downlink from LI; however, a spacecraft at Venus incurs a factor of ~900-30,000 drop in 
data rate relative to LI. Methods to radically increase data rates such as optical spacecraft 
communications systems have yet to be adopted for routine use by flight projects. 

The signihcant decrease in data rate from Venus means that full-frame images cannot 
be downlinked without lossy compression, increasing the difficulty of discovering NEOs at or 
near the sensitivity limit. While the “windowing” method was used on the Kepler mission 
to select only small groups of pixels around targets of interest, this technique requires a 
priori knowledge of the targets’ positions. However, when searching for new NEOs, target 
locations are not known. “Windowing” to reduce data volume for an NEO survey therefore 
requires performing all source extractions on board the spacecraft. Given limitations on 
memory and processor capacity, the ability to compute backgrounds and identify transient 
instrumental artifacts is greatly reduced. Such techniques have not been demonstrated to 
produce calibrated sources at low SNR. Moreover, survey projects such as the WISE and 


Two-Micron All-Sky Survey projects (Skrutskie et ah 2006) benefit substantially from the 


ability to reprocess science data using optimally estimated calibration products derived after 


the survey data were collected (Cutri et ah 2012). If full-frame images are not collected, then 


reprocessing to apply best calibration products cannot be carried out, leaving many potential 
detections lost and ruling out the possibility of improved photometry and astrometry. 

A Venus-trailing IR telescope, at 0.7 AU from the Sun, is further stressed by the in¬ 
creased heat load on the thermal system, requiring active cooling for a cryogenic instrument 
whose focal planes must operate at ~35 K. The overall design and system complexity of a 
spacecraft that must now operate far from Earth rather than remaining relatively nearby is 
increased. 

To compare the asteroid detection capabilities of surveys operated from LI and Venus- 
trailing orbit, we have carried out detailed simulations of these two scenarios. Two simulation 
approaches were taken. First, we performed a detailed simulation of a survey observing a 
limited region of sky over two years. This approach includes the creation of spatially static 
sources (primarily stars and galaxies) based on extrapolation of WISE and Spitzer data to 
NEOCam’s predicted sensitivity limits as well as modification of the WISE Moving Object 


Processing System (WMOPS; Cutri et ah 2012 Mainzer et al. 2011a) to accommodate the 
cadence required to discover NEOs without the need for ground-based follow-up (i.e. self- 
follow-up). The simulated near-Earth asteroid (NEA; near-Earth comets were not included 


in this study) population is based upon the latest estimates from NEOWISE (Mainzer et al. 
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2011a, 2012b), with orbital elements taken from Bottke et al.[(2002) as implemented by Grav 


et al. (2011a). The second approach did not include static sources, but instead simulated 
a survey that covers the entire viewable sky over a six-year period. Taken together, 
these two approaches allow new insights to be gained into the potential capabilities and 
optimization of a next-generation IR survey. 


Both simulations follow this outline: 1) A frame-by-frame list of pointings for the survey 
is generated; 2) mission parameters are input (e.g. instrument FOV including detector gaps, 
sensitivity, slew time, downlink time, etc.); 3) a synthetic asteroid population that models 
expected numbers, orbital elements, and physical properties (e.g. diameter, albedo, etc.) 
is generated; 4) the predicted positions, infrared fluxes, and visible fluxes of each object in 
the population model are computed for each exposure in the simulated survey over six years 
using two thermal models; 5) objects that are detected the requisite number of times at the 
appropriate cadence are tallied. The population model is generated by randomly drawing 


orbital elements from the 

Bottke et al. 

(2002 

) model and assigning diameters and albedos 

randomly from the distributions given in 

Mainzer et al. 

(2012b 

) for Atens, Apollos, and 


Amors. The crucial difference between the limited-region and all-sky surveys (other than 
the area covered) is that in the limited-region simulation uses the modihed WMOPS system 
to link detections together, allowing its efficacy to be evaluated. 


2. Survey Simulation Assumptions 

In this section, we describe the assumptions that were made for both the limited region 
and all-sky survey simulations. 

We assume that both LI and Venus-trailing missions employ similar instruments with 
identical performance characteristics (i.e. telescope aperture size, detector performance, slew 
times, etc.). Both instruments are assumed to cover the sky at the same rate, and both must 
obey the same restrictions on their cadence: that is, while each instrument has different 
viewing constraints owing to their different orbits, the sequence of observations in both 
position and time required to successfully discover NEOs is identical. Even though there 
is a factor of ~900-30,000 decrease in data rate for the Venus-trailing mission compared to 
LI (necessitating complex on-board data processing and destructive data compression that 
will degrade performance), we nevertheless assume that data can be extracted and processed 
to the same sensitivity limits and with the same astrometric and photometric precision as 
the LI mission. No penalty for lossy data compression was assumed. If sources cannot be 
reliably extracted down to low SNR, then performance will be degraded accordingly. 




















2.1. Instrument Performance Characteristics 


Instrument: The requirements for both missions are summarized in Table 1. Both 
missions (LI and Venus-trailing) employ a 0.5 m telescope with a 1.85° x 7.77° held of view. 
Both telescopes are assumed to be diffraction-limited in their bands, and both require 30 sec 
to slew between helds and settle sufficiently for precision pointing. Total integration time 
at each pointing is taken to be 180 sec in an effort to maximize sensitivity and minimize 
inefficiency lost to slew time. Both missions employ a 2 x 8 mosaic of 1024 x 1024 10 
fim HgCdTe detectors bonded to Teledyne Imaging System’s HAWAII IRG readout circuits 
similar to those described in McMurtry et ah (2013), as well as a 1 x 4 mosaic of 5 fim 
HgCdTe detectors in a 2048^ format. The spacing between the arrays is shown in Figure 
Both LI and Venus-trailing missions are assumed to have two bandpasses that simultaneously 
image the same region of the sky using a beamsplitter: 4 - 5.2 fim (referred to as NCI) and 
6 - 10 /rm (NC2). 


CHANNEL NCI: 4.0-5.2 pm 


~3mm ~6 mm 

,,, CHANNEL NC2: 6-10 urn 

IS 4S7 mm ^ 

£ 

£ 


'“S mm 

Fig. 1.— Schematic showing the layout of detector arrays in mosaics for the two channels 
(4 - 5.2 pm, NCI, and 6-10 pm, NC2). The NCI channel consists of a mosaic of 2048 x 
2048 HgCdTe arrays, while the NC2 arrays are 1024 x 1024 format. Both channels image 
the same held simultaneously using a beamsplitter. The gaps between detector arrays were 
included in both all-sky and test region simulations. 





The optical systems for both LI and Venus-trailing missions are taken to be diffraction- 
limited at 4 pm. The point spread function (PSF) full width at half-maximum (FWHM) 
for both missions is ~4 arcsec in channel NC2, and the plate scale is set to 3 arcsec/pixel 
to maximize areal coverage at the expense of some loss of sampling. The PSF will be 
undersampled for the NCI channel. Astrometric uncertainty is taken to be ~0.5 arcsec l-cr. 














- 9 - 


Table 1. Requirements for both an LI and Venus-trailing survey. Values for NCI and 
NC2 channels are separated by Read noise is specihed for correlated double sampling 
(CDS). *Since point source sensitivity varies with ecliptic latitude, longitude, and 
heliocentric distance, the point source sensitivity is specihed near the midpoint of the 
viewing zones for each survey: (A, (3) = (0°, 90°) for LI, and (0°, 180°) at 0.7 AU 
heliocentric distance for the Venus-trailing survey. 


Requirement 

LI 

Venus-trailing 

Telescope aperture (m) 

0.5 

0.5 

Wavelengths (pm) 

4 - 5.2; 6 - 10 

4 - 5.2; 6 - 10 

Field of view (°) 

1.85 X 7.77 

1.85 X 7.77 

Slew time (s) 

30 

30 

Dwell time (s) 

180 

180 

Minimum background current (electrons/s) 

~300 

~300 

Detector format 

1x4 2048^ mosaic; 

1x4 2048^ mosaic; 


2x8 1024^ mosaic 

2x8 1024^ mosaic 

Detector type 

HgCdTe 

HgCdTe 

Pixel size (pm) 

18; 18 

18; 18 

Image FWHM (arcsec) 

3; 4 

3; 4 

Dark current (electrons/s) 

<5; <200 

<5; <200 

Read noise (CDS) 

15; <30 

15; <30 

Quantum efficiency (%) 

>60; >60 

>60; >60 

Point Source Sensitivity* (pJy, 5 -(t) 

50; 150 

50; 150 

Viewing zones (solar elongation) 

±(45° -125°) 

180±75° 

Viewing zones (ecliptic latitude) 

±41.9° 

±41.9° 
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Zodiacal Background: In the 6-10 /im band, an instrument’s optical system and 
detectors must be cooled to achieve natural background-limited performance. The dominant 
source of natural background is thermal emission from zodiacal dust. In the 4 - 5.2 /xm 
band, for an adequately cooled instrument, sensitivity is limited by a combination of dark 
current, read noise, and zodiacal background flux. The zodiacal background is represented 
by a model that varies with ecliptic latitude and longitude based on the fluxes tabulated by 


Wright (1998) and Gorjian et ah (2000). As described in Leinert et ah (1998), the density 


of zodiacal dust increases as 1/r, where r is heliocentric distance. The requirements for 
detector dark current and read noise derive from the minimum zodiacal flux. 


2.1.1. Orbits and Viewing Constraints 


The asteroids likely to be the most hazardous and require the lowest amount of energy 
to reach from Earth (i.e. those with the most Earth-like orbits) are most readily found 
by searching the space around the Earth’s orbit (Chesley & Spahr 2004). Therefore, it is 
desirable to design an optical system for a telescope located at the Earth-Sun LI point to 
be able to look as close to the Sun as possible in order to subtend the largest fraction of 
Earth’s orbit (Figure]^. Similarly, the Venus-trailing mission must be designed to maximize 
the instantaneous viewing area. 



Fig. 2.— Viewing constraints for both the Earth-Sun LI Lagrange point (left) and the 
Venus-trailing (right) surveys as shaded areas. The orbits and locations of Mercury, Venus, 
Earth and Mars are shown. 
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The LI mission is shielded from the Sun sufficiently to allow it to point anywhere from 
45 to 125° solar elongation in either direction. IRAS could observe as close as 60° from the 
Sun ( Beichman et al.||1988 ), and the Hubble Space Telescope has carried out observations as 
close as ~ 45° solar elongation (Na & Esposito 1995). The Venus-trailing mission’s viewing 
zone is centered at opposition and can point up to 75° away in elongation. Despite its larger 
instantaneous viewing zone area, the Venus-trailing mission cannot cover more “fresh” sky 
per unit time than the LI survey, and it must execute a survey cadence that supports NEO 
discovery. 


2.2. Survey Cadence 


The Minor Planet Center (MPC), a division of the International Astronomical Union, 
arbitrates discoveries of new minor planet candidates. In order for an object to be declared 
“discovered,” it must be observed over several days. While this is sufficient to classify most 
objects as NEOs or more distant orbits, this observational arc is generally insufficient to 
allow objects to be recovered by observers at future apparitions using conventional wide- 
held imagers; such objects are often effectively lost. For an object to be recoverable at 
its next apparition, it is generally necessary to extend its observational arc out to at least 
10-25 days (preferably 60 to 90 days), depending on the observatory’s astrometric accuracy. 
That is, astrometric uncertainty can be no more than a few degrees off its prediction position 
after ~3-10 years following its hrst detection. It should be noted that observational arcs of a 
few weeks in duration are generally insufficient for formal hazard evaluation, and in general 
observational arcs spanning more than one opposition are required to reduce the orbital 
uncertainty sufficiently to allow location of the object within a few arc seconds. However, 
even the shortest arcs described here, spanning ~10 days, will allow linkage at subsequent 
return visits. In essence, this is routinely accomplished by the MPC today as demonstrated 
by past recoveries and linkages of (719) Albert, 1937 UB (Hermes), and even 1954 XA. 


Unlike ground-based surveys, space telescopes at either LI or Venus-trailing orbits can 
spend much of their time surveying the region of sky that is in the daytime sky for ground- 
based observers. Therefore, the ability to perform “self follow up” is essential because ground- 
based observers cannot be relied upon regularly for the short-term follow up required to 
determine orbits securely. These space-based surveys are designed to reach dimmer limiting 
magnitudes than most existing ground-based follow-up stations (see Section 4.1 and Figure 
22). Future surveys must employ a cadence that ensures detection, follow-up, and accurate 


orbit determination are built into routine operations. 


Experience at the MPC has shown that a >20 day arc with approximately a dozen or 
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more observations spaced over time is required to fit an orbit with sufficiently low astrometric 
uncertainty to allow next-apparition recovery. The cadence we tested for both the LI and 
Venus-trailing surveys is illustrated in Figures 3-6 (spherical geometry has been neglected 
for simplicity). The sequence begins with a set of four images taken ~1-1.5 hours apart (a 
“quad”), spanning ~9 hours (Figure [^. This cadence has been demonstrated by ground- 
based surveys such as Catalina to be robust against false linkages. The survey starts at 
the maximum solar elongation and lowest ecliptic longitude. Each step is taken along the 
shorter direction of the detector mosaic to minimize slew time. Next, this loop pattern is 
repeated an additional hve times, each at increased ecliptic latitude up to the maximum, 
over the course of ~2.25 days (Figure]^. For the LI survey, this block of helds is surveyed 
over ~6 days on one side of the Sun (Figure [^. Next, an identical pattern is surveyed for 
the following ~6 days on the other side of the Sun, after which point, the survey returns 
to the original position and repeats the sequence a total of 11 days later (Figure]^. Both 
the LI and Venus-trailing surveys cover sky at the same rate, and both result in ~22 day 
observational arcs for a large fraction of the NEOs observed. 




1 Loop = 

126 min; 
repeat 4 
times to get 
4 detections 
in ~9 hours 


= Starting field 


Fig. 3.— The basic element of the survey cadence is a “quad,” which consists of four 
observations taken over ~9 hours. The green rectangle indicates the starting position. 
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+41.9° latitude 


-41.9° latitude 



4 loops, 
''44-53 hrs 


4 loops, 
■''36-44 hrs 


4 loops, 
-27-36 hrs 


4 loops, 
-18-27 hrs 


4 loops, 
-9-18 hrs 


4 loops, 
-0-9 hrs 


2.25 days 


A 


J 


Fig. 4.— The loop pattern is repeated, increasing in ecliptic latitude, over ~2.25 days. The 
green rectangle indicates the starting position. 
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~6 days 


Fig. 5.— The survey pattern over ~6 days for the LI survey. Green indicates the starting 
held. 


Total of ~12 days from start 
to finish on each area 


We shift over to account for the 
Sun's motion on the sky and 
begin the entire pattern again. 
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(110°, -41.9° 


(45°, -41.9°) 


(45°, 41.9°) 
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(45°, -41.9° 


(110°, -41.9°) 


Fig. 6.— The entire cadence takes ~22 days to execute, switching back and forth on either 
side of the Sun for the LI survey, and proceeding in a continuous sweep for the Venus-trailing 
survey. 
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2.3. Synthetic Asteroid Population Models 


For both all-sky and test region simnlations, we considered NBAs as small as 140 m 
in effective spherical diameter. The popnlation models inclnded Atens (NBAs with aphelia 
Q >0.983 AU and semi-major axes a <1.0 AU), Apollos (NBAs with a >1.0 AU and perihelia 
q <1.017 AU), Amors (NBAs with a >1.0 and q >1.017 AU), and interior-to-Barth objects 
(IBOs, or Atiras; NBAs with a <1.0 AU and Q <0.983 AU). For the all-sky simnlation, the 


synthetic solar system orbital element models of both Grav et ah (2011a) and Greenstreet 


& Gladman (2013) were nsed to generate 25 randomly drawn popnlations of Atens, Apollos, 


and Amors. By rnnning 25 Monte Garlo trials, we average over variations in orbital elements 
and physical properties in randomly generated popnlations. Greenstreet et al. (2012) was 


nsed to generate the IBO popnlations’ orbital elements, as these were not inclnded in the 


Grav et al. 

(2011a 

) model, which is based on the orbital element distribntion of 

Bottke et al. 


(2002). The nnmbers of Atens, Apollos, and Amors for the 25 synthetic popnlations were 


chosen according to the predicted total nnmbers fonnd by Mainzer et al. (2011b, 2012a), for 
a total of ~13,200±1900 NBAs >140 m in each synthetic popnlation. For the test region 
simnlation, only a single simnlated popnlation was generated, containing 12,700 NBAs >140 
m. The simnlation does not at present consider near-Barth comets or long-period comets 
that enter near-Barth space. 


Minimnm orbital intersection distances (MOIDs; Bowell & Mninonen 1993) were com- 


pnted for each synthetic object nsing the same methods that are employed by the JPL 


Horizons system (Ostro & Giorgini 2004). MOIDs were compnted so that the synthetic 


NBOs conld be separated into those that wonld be deemed potentially hazardons asteroids 
(PHAs), depending on their size. PHAs are formally defined as having Barth MOIDs <0.05 


AU and having absolute H < 22 mag, althongh Mainzer et al. (2012a) snggest a diameter- 


based definition, since IR snrveys do not directly sample H. 

The physical properties for each object were randomly assigned according to the distri- 


bntions given in Mainzer et al. (2012a) for size, geometric albedo {pv)^ beaming parameter 


7] employed by the Near-Barth Asteroid Model (NBATM; Harris [r998 ) and infrared albedo 
{piR, the reflectivity at ~3-5 pm) for the Atens, Apollos, and Amors, respectively. Mainzer 


et al. (2012a) fonnd that the Atens, Apollos, and Amors each have different size and albedo 


distribntions. While the definitions of these three gronps are somewhat arbitrary, the differ¬ 
ences in physical properties probably reflect the differing sonrce regions’ properties for each 
gronp. Since only 14 IBOs have been discovered to date, little is known abont their total 


nnmbers or physical properties. Popnlations of IBOs were generated nsing the Greenstreet 


& Gladman (2013) orbital element model. Lacking information on the estimated total nnm¬ 


bers of IBOs and their physical properties, the total nnmbers and size freqnency distribntion 












































16 


were assumed to be similar to the Atens; their albedo distribution was modeled based on 
the Apollos. Figure shows the median size and albedo distributions for all synthetic NEA 
subpopulations. 


We created 25 synthetic populations to lessen the chance that one individual randomly- 
generated set of asteroids could represent an outlier in terms of numbers of objects, orbital 
elements, or physical properties. While the actual numbers, orbital element distributions, 
and size distributions are intended to be measured as a scientihc goal of a next-generation 
survey, this set of model populations allows comparison of two potential survey implemen¬ 
tations (LI and Venus-trailing). 


Fluxes were computed at NCI and NC2 wavelengths for all objects at each time step in 


the simulated surveys using both the NEATM and the Fast Rotating Model (FRM; Lebofsky 
et ah 1978|[Veeder et ah 1989 Lebofsky & Spencer 1989). Objects were modeled as faceted 


spheres (c.f. Kaasalainen et ah 2004) as described in Mainzer et al. ( 2011b|[c ). The NEATM 
assumes a temperature distribution that decreases from the subsolar point to essentially zero 
at the terminator, whereas the FRM assumes a temperature distribution that is uniformly 
distributed in longitude. While the NEATM assumes that the nightside of each asteroid 


contributes no flux, thermophysical models (e.g. Groussin et al. 2011 Delbo & Tanga 2009 


Lebofsky & Spencer 1989) have shown that the assumption of zero nightside flux is not 


always good, particularly for observations at high phase angles where more of the nightside 
is seen. Objects with high thermal conductivity will tend to have more uniform temperature 
distributions. The NEATM also requires that the subsolar point is the hottest surface 
point. Thermophysical models also show that with hnite rotation and thermal inertia, the 
peak temperature is lowered and displaced from the subsolar point, altering the day side 
temperature distribution as well. Small NEOs are known to have a higher fraction of fast 


rotators than larger objects (Warner et al. 2009 Pravec et al. 2008 Pravec & Harris 2007 


2000). It is therefore useful to compare NEATM and FRM results as reasonable bounding 


cases. The survey simulation results for the all-sky LI and Venus-trailing surveys using both 
the NEATM and FRM are given in Section 4.1. 
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Fig. 7.— Mean cumulative size frequency distributions (left) and albedo distributions (right) 
for the 25 synthetic populations of lEOs, Atens, Apollos, Amors, and all NBAs used in the 
all-sky survey simulation. 


3. Test Region Simulations 

In this section, we describe the results of simulations that were carried out for both 
Venus-trailing and LI surveys on a limited test region spanning ~260 square degrees popu¬ 
lated with simulated NBAs, the full catalog of ~600,000 known asteroids and comets, and 
static sources such as stars and galaxies. The purpose of these detailed simulations was to 
test the efhcacy of our moving object detection algorithms in the presence of potentially 
confusing transient noise and static sources. The ~260 square degree held is centered on 
ecliptic longitude/lattitude (-1-140°, latitude 0.0°), bounded by 129 to 151° longitude and -6 
to -|-6° latitude. 

A list of held center positions and spacecraft ephemerides for the hrst two years of 
the Venus-trailing and LI surveys were used as the basis for source list generation using 
the observing cadence described in Section 2.2. A background “static sky” was generated 
for each frame based on source count measurements from previous IR surveys as described 
below. Fluxes for a population of minor planets, including a synthetic NEA population 
and the full catalog of previously known Main Belt asteroids, Jovian Trojans, planets, etc., 
were computed for each frame and added to its source list. Position-time sets for each 
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candidate minor planet detection (“tracklets”) were linked nsing a version of WMOPS that 
was adapted for the specihc needs of the new cadence. This version of the pipeline, dnbbed 
the Experimental Moving Object Processing System (XMOPS), allowed us to estimate the 
tracklets’ reliability and completeness as well as the frequency of cross-linkages between 
inertial sources and other asteroids. The final results for the test region simulations are 
described in Section 3.4. 


3.1. Static Sky 


The static sky consists of astrophysical sources such as stars and galaxies that could 
potentially be confused with minor planets. For our test region simulation, we constructed an 
artihcial “sky” using source catalogs from WISE and the Spitzer Space Telescope’s Infrared 
Array Camera (IRAC; Fazio et ah 2004), along with knowledge of the mid-IR sky (e.g. 


Jarrett et al. 2011| . The survey simulation bandpass of NCI is approximately the same 
as the WISE 4.6 /rm (hereafter referred to as W2) and IRAC 4.5 /rm (IRAC-2) channels, 
whereas the longer band, NC2, has a response that is closer to the IRAC 8 fim (IRAC-4) than 
the WISE 12 fim channel (W3). Hence, we assume that IRAC-4 may be used to construct 
the NC2 sky. However, since IRAC-4 images are only available over limited regions of the 


sky, such as the Spitzer Wide-area Extragalactic Survey (SWIRE; Lonsdale et al. 2003), we 


use the all-sky W2 and W3 densities to predict the IRAC-4 and hence NC2 sky. 

We did a series of tests to determine the best method for extrapolating the static sky 
in the NCI and NC2 bandpasses. We used the SWIRE and WISE catalogs of the European 


Large Area Infrared Space Observatory Survey North region (ELAIS-Nl; Rowan-Robinson 


et al.|[l9^ to explore the relationship between IRAC-2 and W2, and IRAC-4 to W2 and 
W3. Figure]^ (left) shows the source counts in W2 and IRAC-2 from the ELAIS-Nl region. 
The plot includes the expected source counts for foreground stars in the Milky Way, as well 
as background galaxies as determined in the near-infrared M band. W2 and IRAC-2 agree 
reasonably well for W2 < 16th mag (68 /iJy), but at fainter magnitudes WISE becomes 
incomplete and is diminished relative to IRAC-2. The figure also shows the extrapolation to 
fainter magnitudes, a necessary step because WISE is not deep enough to reach the predicted 
sensitivity level in the NCI band. Although the fit to the WISE counts appears reasonable, 
the extrapolated source counts slightly underestimate the IRAC-2 counts at fainter flux levels 
(IRAC-2 appears to be incomplete at about ~20 /iJy). The discrepancy between the W2 and 
IRAC-2 at the faint levels likely arises from both incompleteness (WISE) and the Eddington 
bias, which causes fluxes to be overestimated at the faint end (IRAC-2). 
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W3 [Vega mag] 


Fig. 8.— Left: Expected source counts for foreground stars in the Milky Way (blue solid 
line) and background galaxies (purple solid line) based on W2 (4.6 /im) and IRAC-2 (4.5 /im) 
source counts from the ELIAS-Nl region; both W2 (black solid line) and IRAC-2 (orange 
dashed line) are very similar to NCI (4.6 pm). At the faint end, the red dashed line shows 
the extrapolation to fainter magnitudes. Right: Expected source counts for stars (blue solid 
line) and source counts from IRAC-4 (8 pm; green dashed line) and W3 (12 pm; black solid 
line). The red and orange points show the predicted total source counts for sources in the 
IRAC-4 band (taken to be a proxy for NC2) based on extrapolation from W2, W3 and 
IRAC-4. 


Figure (right) shows the resultant W3 and IRAC-4 source counts for the ELAIS-Nl 
held. The NC2 band (6 - 10 pm) is not closely aligned with that of W3 (which spans ~7.5 
- 16.5 pm), but it is reasonably close to that of IRAC-4 (8 pm). Since we do not have 
IRAC-4 images for most of the sky, we must use the WISE W3 all-sky data to reconstruct 
the IRAC-4 sky. We must also account for the source counts in W2 since the static sky is 
a catalog of bandmerged NCI and NC2 sources. We therefore use W2 to predict the W3 
source counts, then use our knowledge of W3 to predict the IRAC-4 counts. The band- 
to-band differences between W2 and W3 are large, so the relative fraction of stars (which 
are blue in color at these wavelengths) and galaxies (which are red) makes a signihcant 
difference; i.e. helds dominated by foreground Milky Way stars will have a different color 
than regions dominated by extragalactic sources. Figure]^ (right) shows that the W3 counts 
can be reasonably reproduced by using W2 as the primary seed, and then the IRAC-4 source 
counts can be predicted from W3 accordingly. Although the experiment on this region of sky 
was successful, more work remains to be done to verify that this method is robust against 
large color differences caused by a different mix of stellar and extragalactic sources. 

Source counts for the test region simulation. Using the methods described above, we 
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generated the static sky for the entire ~260 sqnare degree test region simnlation. Fignre|^ 
(left) shows the resnltant W2 connts. The W2 connts become incomplete for W2 >16 mag 
(68 /iJy); htting to the faint bins, the resnlting £t may be nsed to reconstrnct the W2 connts 
for all flux levels, notably those that are near the NCI sensitivity limit. To simulate the NC2 
static sky, we use W2 and W3 to extrapolate the source counts for IRAC-4, which serves as 
a proxy for NC2 as shown in Figure]^ (right). The reconstructed W3 population agrees well 
with the observed W3 source counts down to the limit at which W3 becomes incomplete. 
This good agreement gives conhdence that the extrapolated IRAC-4/NC2 source counts 
represent the NC2 static sky accurately. 

As a hnal step in generating the static sky catalog for the test region, random right 
ascension and declination positions were assigned to “new” sources that are fainter than the 
real W2 limit. Position uncertainties were modeled as a function of intensity. 



20 15 10 5 

W2 [Vega mag] 



W3 [Vega mag] 


Fig. 9.— Left: A £t to the W2 source counts (blue dashed and black solid lines, respectively) 
produces the extrapolated NCI source counts (orange dashed line). The W2 counts become 
incomplete at ~16th magnitude. Right: W2 and W3 are used to reconstruct the IRAC-4 
counts (dashed magenta line), which are used as a proxy for the NC2 sky (see Figure]^ for 
method details). 


3.2. Moving Objects 

In order to create a population of non-NEAs that could create potentially confusing 
cross-links to NEA candidates, known asteroids were added to the XMOPS test frames. The 
entire catalog of known asteroids was downloaded from the MPC, and objects with perihelia 
>1.3 AU, semi-major axes <34.19952, and eccentricity <1 were injected into the synthetic 
source lists. A single synthetic NEA population consisting of 12,700 objects larger than 140 
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m was generated as described in Section 2.3 above. 


3.3. Experimental Moving Object Processing System 


Sonrce lists comprised of static sky backgronnds and synthetic MBA and NEAs were 
generated for each frame in the snrvey simnlation that fell within the test region. A two year 
interval was considered, from Jnne 1, 2012 to Jnne 1, 2014. Asteroid flnxes in the NCI and 
NC2 channels were generated using the NEATM. Simulated position and magnitude were 
errors added to each source. 


WMOPS was modified slightly to operate on source lists generated according to the new 
survey cadence (three sets of four quads spanning ~22 days). The resulting system, XMOPS, 
was run on the two year’s worth of simulated source extractions and compared to the known 
number of NEOs that appeared in each frame to compute an estimate of the efficiency with 
which tracklets could be found. XMOPS provides the same core functionality as WMOPS, 
with the major differences being in the way in which regions are selected and grouped for 
processing. The test region of sky was broken into smaller 3x3° subregions (“patches”) that 
overlapped by ~50% or more. Stepping along in 1.5° increments, each patch was checked to 
see if it had been imaged at least four times over a 24 hour interval. If there were at least 
four coverages in 24 hours, the WMOPS routine that identifies inertially fixed sources with 
SNR>4 was run to remove sources that repeated at the same location from the source lists. 
As described in Section 2.2, a minimum of four detections spanning ~9 hours is required 
to form a “quad”. A later threshold at SNR=6 was applied to the tracklet-detection input 
lists, and the tracklet efficiencies were used to inform the lower SNR samples (see Section 
2.5.1). 


The FindTracklets routine (|Kubica et al. 2007| was run to link pairs of candidate tran¬ 


sient sources, similar to the methods employed by WMOPS (Mainzer et al. 2011a Cutri 


et al. 2012). The Kubica et al. (2007) method uses hierarchical data structures called k-d 


trees to recursively partition the sources that could potentially be linked into smaller sub¬ 
sets, reducing the search time proportional to plogp (where p is the source sky-plane density) 
instead of p^. The CollapseTracklets algorithm (Myers et al. 2008) was used to link the de¬ 
tection pairs into candidate lists of moving object tracklets. Minor adjustments were made 
to the parameters limiting how closely spaced in time the detections had to be as well as the 
maximum allowable velocity (Figure 10). Since WMOPS relied upon an “eyes-on” human 
quality assurance step in which new candidate object tracklets were visually inspected, but 
a similar system was not implemented in this simulation, the maximum allowable velocity 
was reduced slightly to increase the reliability of tracklet linkages. An example of candidate 
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tracklets generated by XMOPS is shown in Fignre The incidence of cross-linkages be¬ 
tween two different asteroids and stationary sources that were not correctly identihed by the 
static sky rejection algorithms were evaluated by comparing XMOPS outputs to the list of 
true tracklets. 


350 




— LI 

— Venus-trailling 



Velocity (deg/day) 


Fig. 10.— The average apparent on-sky velocity of NBAs in two years’ worth of frames 
simulated in the ~260 square degree test region searched with XMOPS. NBAs detected by 
the Venus-trailing survey tend to have lower on-sky angular velocities (red line) because they 
are observed at greater heliocentric distances (this also makes them fainter). Detections can 
only be linked if the apparent velocity is low enough to avoid signihcant trailing. 
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Fig. 11.— A set of example candidate tracklets with four or more detections generated for 
one patch by XMOPS (red arrows indicate the start and end points of the tracklets). Empty 
blue “o” marks show the simulated detections that match the position, time, magnitude, 
and coverage criteria and were not identihed as stationary sources. Black “x” marks show 
the expected, true known solar system objects that match the time/position window as well 
as the magnitude and coverage thresholds. An “x” without an “o” is a true, detectable 
solar system object that did not survive stationary object rejection, and an “o” without an 
“x” is a stationary source that survived stationary object rejection. Any “x” or “o” that 
is not covered by a red arrow was not linked into a tracklet. All but two of the red arrows 
represent Main Belt asteroids; the sole NEA in the frame is highlighted by the yellow circle. 
An example of an erroneously cross-linked tracklet is shown by the nearly vertical red arrow 
on the extreme left. 
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3.4. Results 

For the LI mission, there were 1,123 NBAs that appeared with the minimum number of 
coverages (>4), but only 341 of those appeared above the SNR=6 cutoff. Of these objects, 
313 were successfully linked into one or more tracklets, yielding a linking success rate of 
92%. The 112 NBAs that were detected in two or more tracklets represented 87.5% of the 
128 NBAs detected above the SNR=3 threshold. These linking success rates are nearly 
identical to the results for the Venus-trailing survey, which linked 160 out of 170 NBOs 
above the SNR=6 cutoff into one or more tracklets (94%), and 47 out of 55 NBAs into two 
or more tracklets (85%). Only 385 NBAs appeared above the coverage threshold for the 
Venus-trailing survey during the two year simulation; the number is considerably decreased 
relative to the LI mission because the observatory is circulating at Venus’ orbital period 
rather than Barth’s, and the XMOPS portion of the simulation does not cover the entire 
sky. 

The nearly identical tracklet linking rates for both surveys suggests that there is no 
particular advantage or disadvantage to either orbit insofar as the ability to link tracklets is 
concerned, assuming that sources can be reliably extracted to low SNR. A limitation of the 
XMOPS simulation at this point is that it did not consider the impact of spurious detections 
of artifacts and noise on the number of possible combinations that can create a tracklet. 
If spurious detections cannot be reliably identihed and excluded from the list of transient 
sources, they will overwhelm the pipeline with false cross-linkages. The ability to correctly 
identify and exclude transient artifacts such as stray light, cosmic rays, and noisy pixels is 
essential. 


4. All-Sky Simulation 


A second set of simulations was created to estimate the theoretical best performance that 
an advanced survey could achieve over the entire viewable sky for a six year survey for both 
LI and Venus-trailing missions. These models did not make use of XMOPS in the presence 
of a static sky. Instead, we produced a list of every held pointing of the six year surveys and 
compared this to the positions of populations of simulated objects. We generated synthetic 
orbital elements for each simulated object, and then numerically integrated the positions of 
these objects forward over the duration of the survey under the gravitational inhuence of 
the Sun, the eight planets, and the Moon with one hour timesteps. We used the SWIFT 
numerical integrator ( ]Levison fc Dunclin 1994) implementing the Bulirsch-Stoer integration 
method on JPL’s high performance computing facilities to calculate the positions of each 
object over the course of our simulation. Positions of each object over time were then 
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compared with the date and coverage of each held pointing, and along with distance from 
the telescope nsed to determine detectability. 


As described above, the natnral sky backgronnd was assnmed to vary according to the 
models and measnrements of Wright (1998) and Gorjian et ah (2000) as a fnnction of ecliptic 
latitnde and longitnde. The Venns-trailing survey’s zodiacal background was assumed to 
increase as a function of heliocentric distance according to Figure 55 in Leinert et ah (1998). 
Both Venns-trailing and LI surveys were assumed to have the properties given in Table 1 and 
were assumed to carry out the cadence described above. The sensitivity of the instrument 
was computed over a grid of ecliptic latitudes and longitudes with 5° steps spanning the 
helds of regard of both Venns-trailing and LI surveys (±41.9° latitude for both surveys; 45 
- 125° solar elongation on either side of the Sun for the LI survey at ~1 AU heliocentric 
distance, and ±75° in elongation centered at opposition at 0.7 AU heliocentric distance for 
the Venns-trailing survey). For each NEA’s position in each pointing during the surveys, the 
sensitivity was determined by interpolation between grid points in the sensitivity model. 


In these simulations, an object was counted as detected if its predicted flux was above 
the SNR threshold and its velocity fell within an allowable range. The upper and lower 
apparent velocity limits were determined by trailing losses and the instrument’s astrometric 
accuracy for faint sources, respectively. This survey simulation did not account for confusion 
with inertially hxed sources such as stars and galaxies, confusion with optical or electronic 
artifacts, or more distant asteroids that could potentially be confused with NEOs. Therefore, 
these simulated surveys should be regarded as predicting the “best-case” of an advanced IR 
survey appropriate to low-confusion regions of the sky. 


The all-sky survey simulations for both LI and Venns-trailing missions were assumed to 
start on June 1, 2014 and end on June 1, 2020. The all-sky survey assumes that sources are 
extracted down to SNR=5. As described in Section 2.3, a total of 25 synthetic populations 
were generated. Using 25 randomly drawn populations ensures that the results are less likely 
to be affected by outliers in orbital elements or physical properties. The all-sky simulations 
were performed using NEA fluxes computed with the NEATM and fast rotating thermal 
models for both Venns-trailing and LI surveys. 


Since an advanced survey will build upon the progress made by existing ground-based 
facilities, it was necessary to develop a method of estimating how many NEAs in a given 
size range would have already been discovered by these projects. A synthetic survey of 
ground-based detections beginning in 1985 and proceeding in hve year increments to the 
assumed end date of the advanced survey (June 1, 2020) was created. Visible light fluxes 
were computed for each synthetic NEA between 1985 and 2014. Objects were declared 
“found” by the visible light surveys if they appeared in the night sky with solar elongations 
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> 120°, ecliptic latitudes < 40°, and apparent magnitudes brighter than specihed limits for 
each survey epoch. The limiting magnitude of each epoch was tuned to the sensitivity of 
each successive generation of NEO survey telescopes from 1985 to the present day. The 
resulting H magnitude distribution compares very well to the known population. The all¬ 
sky simulations assume that the existing suite of ground-based surveys (e.g. Catalina and 
PanSTARRS) continue to operate throughout the duration of the LI and Venus-trailing 
IR missions. The comparison between a space-based IR mission and an enhanced ground- 
based survey will be the subject of future work. Ground-based and space-based missions 
complement each other, since ground-based telescopes are particularly sensitive to Amors, 
whereas the space-based surveys detect a higher fraction of lEOs, Atens, and Apollos as 
described in Section 4.1 below. While diameters can be determined using IR fluxes alone 
once a reliable orbit has been established, visible albedos can only be determined if objects 
are detected by both IR and visible light surveys, signihcantly enhancing the value of both 
datasets. 


4.1. Results 

Figures and show the year-by-year integral survey completeness results for 140 m 
objects for NEATM and FRM, respectively. Integral survey completeness is the total fraction 
of all objects greater than or equal to a certain size that have been detected. Estimates of 
the number of objects that were observed by the space-based surveys but were previously 
discovered by the ground-based surveys are included in the total. The model of projected 
optical survey performance is a simple one and needs rehnement, but it is adequate for the 
purpose of comparing the LI and Venus-trailing surveys. The top rows in both hgures show 
the integral survey completeness for NEAs detected four times over ~9 hours (an individual 
tracklet); the bottom rows show the integral survey completeness for NEAs detected with 
>22-day observational arcs. At present, our simulations include only a relatively simplistic 
method for linking detections into tracklets across 22-day intervals in an attempt to replicate 
the function of the MFC. We were not able to simulate more complex linkages, e.g., a ten-day 
observational arc followed by another set of ten observations carried out six months later. 
Such a set of hypothetical observations is very likely to be linked by the MFC, resulting in 
a discovered object with a well-determined orbit. Therefore, the >22-day arc results likely 
represent a conservative view of the surveys’ capabilities. The objects with four observations 
spanning ~9 hours represent an upper limit, since objects are typically not designated as 
discovered until after they have received observations on two nights separated by >12 hours. 
The one-day and 22-day results therefore represent bounding cases. 
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Fig. 12.— The integral survey completeness for populations of NBAs larger than 140 m 
in diameter versus time, with fluxes computed using the NEATM. Solid lines represent the 
LI survey; dashed lines represent the Venus-trailing survey. Dashed gray lines represent 
2/3 and 90% integral completeness limits, respectively. Left plots: all NBAs (black lines) 
and PHAs (red lines), respectively. Right plots: lEOs (cyan lines), Atens (magenta lines), 
Apollos (green lines), and Amors (blue lines). 
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Fig. 13.— The integral survey completeness for populations of NBAs larger than 140 m 
in diameter versus time, with fluxes computed using the FRM. Color coding is the same as 
Figure [T^ 


Figures and show that the two surveys achieve similar integral survey completeness 
levels for NBAs >140 m. The LI survey outperforms the Venus-trailing survey for Atens, 
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Apollos, and lEOs. The Venus-trailing survey discovers more Amors because it observes 
around opposition. However, Amors are significantly less likely to make close Earth ap¬ 
proaches than Atens and Apollos, decreasing their potential hazard as well as their propensity 
to make suitable rendezvous targets. The LI survey slightly outperforms the Venus-trailing 
survey for PHAs in this size range. The error bars in Figures and [T^ were generated from 
the standard deviation among the 25 simulated populations. The rather small error bars 
indicate that the estimates presented are not significantly affected by sampling error (a fact 
that could not have been reliably known without having done more than one simulation). 

Figure [T4| shows the NEA differential completeness for both surveys using the FRM after 
6 years. Differential completeness does not depend on an assumed size-frequency distribution. 



Diameter (km) 


Fig. 14.— The differential survey completeness for populations of NEAs (black lines) 
and PHAs (red lines) larger than 140 m in diameter for a six year survey, with fluxes 
computed using the FRM. Solid and dashed lines represent the LI and Venus-trailing surveys, 
respectively. 
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Figures [T^ and [T6| show the orbital element distributions for one of the 25 input synthetic 
populations of lEOs, Atens, Apollos, and Amors, along with the orbital element distributions 
of the objects detected with >22 day observational arcs after five years for the LI and Venus 
trailing surveys. As expected based on the results of Figures 12 and 1^, the LI survey is 
more efficient at detecting lEOs, Atens, and Apollos, and the Venus-trailing survey is better 
at detecting Amors. 



Fig. 15.— The distribution of semi-major axis vs. eccentricity for one of the 25 input 
synthetic NEA subpopulations (top row), compared with the distributions of objects detected 
with >22 day observational arcs that were found by the LI survey (middle row) and the 
Venus-trailing survey (bottom row) after Eve years. 
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Fig. 16.— The distribution of semi-major axis vs. inclination for one of the 25 input 
synthetic NEA subpopulations (top row), compared with the distributions of objects detected 
with >22 day observational arcs that were found by the LI survey (middle row) and the 
Venus-trailing survey (bottom row) after Eve years. 


Figure 17 shows the average number of detections per object over the course of a six- 
year survey for both LI and Venus-trailing surveys. Many NEAs receive approximately two 
dozen or more detections, with a large fraction receiving >50 collected from multiple viewing 
geometries. Figures [T^ - show sample lightcurves and viewing geometries for two objects, 
one with an average number of observations, and one with >50 observations. Most objects 
are seen at a minimum of two different viewing geometries, which is useful for thermophysical 
modeling. 
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Fig. 17.— The number of detections per NEA for LI (black lines; heavy and lighter weight 
lines represent NEATM and FRM results) and Venus-trailing surveys (red lines) after hve 


years. 










NC2 Flu 


173.1 
173.0 
172.9 

Fig. 18.— Sample observations of a synthetic Aten detected at two epochs spanning ~24 
days (black diamonds); this cases represents a typical object detected by the LI survey. 
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Fig. 19.— Viewing geometry of the epochs at which the Aten shown in Figure was 
detected. The orbits of Mercury, Venus, and Earth are shown for reference; the positions 
of Earth and the asteroid are shown as cyan and magenta points, respectively. The orbit 
of the Aten is shown in green; magenta lines represent the line-of-sight between the Earth 
and the asteroid at each epoch. This object has a 0.58 year orbital period, so it is observed 
by the LI survey a nearly the same point in its orbit ~423 days after its first apparition as 
indicated by the dates shown in the figure, albeit from a different angle. 
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Fig. 20.— Sample observations of a synthetic Apollo detected 44 times over ~6 months 
(black circles; each circle contains ~3-13 individual observations). The cyan points show 
the flux of the object when it passes through the LI survey’s held of view but is not bright 
enough to be detected. 



36 



AU 


Fig. 21.— Viewing geometry of the two epochs at which the Apollo shown in Figure was 
detected. This object has a 2.3 year orbital period. Color-coding is the same as Figure [T^ 


The distribution of visible magnitudes for the LI survey is shown in Figure both at 
the epoch of the first observation of each object and averaged over each tracklet throughout 
the hrst year of the survey. Visible magnitudes are computed from absolute magnitude H 
using the geometric visible albedo and diameter assigned to each synthetic object and the 
relationship 


H = -5- logio 


1329 


( 1 ) 
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where D is the diameter in kilometers (Bowell et ah 1989). The average visible magnitude 
of objects detected by the LI survey is V~25 mag. 


Figure [23] shows the absolute values of the synodic periods of NBAs detected by both 
LI and Venus-trailing surveys, where synodic period S measured in years with respect to 
Earth is defined as 


Figure shows that the Venus-trailing survey achieves essentially identical performance in 
terms of the ability to detect NBAs of various synodic periods. At the integral completeness 
levels achieved, the survey simulations hnd no “hidden” population of NBAs with decades- 
long synodic periods that are unobservable to the LI survey, since the NBAs librate into the 
helds of regard of the LI survey and do not stay directly behind the Sun. 


The net result is that, even assuming there are no penalties to sensitivity or reliability 
from the lossy data compression needed from Venus-trailing orbit, the LI survey outperforms 
the Venus-trailing survey slightly in terms of integral survey completeness and number of 
detections per object for NBAs larger than 140 m. 
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Visible Magnitude 

Fig. 22.— The distribution of visible magnitudes for the LI survey. The brightness at the 
time of the hrst observation for each object is shown (green line), along with the average 
brightness for each object’s tracklets during the hrst year of the survey (black line). 
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Synodic Period (years) 


Fig. 23.— The distribution of the absolute values of synodic periods for all NBAs in one 
of the 25 synthetic populations is shown (black lines) for both LI (top) and Venus-trailing 
surveys (bottom). The red lines indicate the synodic period distribution for all objects found 
with >22 day observational arcs by each survey after hve years. 


We also studied the case of varying the range of solar elongations for the LI survey, 
including scanning through opposition (which requires moving the spacecraft to L2 to avoid 
viewing the Earth). Assuming that such a survey covers sky at the same rate and must 
carry out the same cadence, surveying through opposition from 110° solar elongation detects 
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only 34% of PHAs after six years and no lEOs. This poor result is not surprising, given 
that Chesley & Spahr (2004) found that PHAs and potential impactors are preferentially 
found at lower solar elongations. Moreover, surveying through opposition presents serious 
engineering obstacles in terms of keeping sunlight off sensitive thermal and optical surfaces. 


There are several areas for improvement in the simulation. First, the synthetic source 
lists created for use with XMOPS did not include models of spurious detections. Experience 
with WMOPS has shown that unidentihed artifacts can severely impact tracklet reliability, 
including those artifacts associated with bright moving objects (e.g. Jupiter and Mars). 
Artifact models based on the HgCdTe arrays used by WISE could be created and injected 
along with models of spurious detections such as diffraction spikes and ghosts. Second, the 
XMOPS portion of the simulation could be extended from its current ~260 square degree 
footprint spanning two years to cover the entire sky over a longer duration. Third, the 
population of co-orbitals such as Earth Trojans and NEAs in so-called “horseshoe” orbits 
must be estimated and included in population models. Finally, it may be possible to develop a 
survey cadence that is better optimized for NEA detection and characterization. The cadence 
described here represents an initial result that could potentially be further optimized. 


4.2. Orbit Determination 


A critical parameter in determining the success of the trial cadence tested for both 
surveys is the ability to compute reliable orbits from the simulated observations. To that end, 
we evaluated the orbits fit to a set of candidate tracklets resulting from both LI and Venus- 
trailing surveys. The determination of the success or failure of any asteroid observational 
cadence is relatively straightforward. Objects must be linked with a high enough reliability 
and observed over an arc of sufficient length that the vast majority of the objects can be 
identihed at subsequent appearances. Orbit fits to a set of candidate tracklets resulting from 
both LI and Venus-trailing surveys using the standard, existing tools employed for real-time 
linking by the MPC were evaluated. 


The method is as follows: First, a Vaisala orbit is computed for each tracklet. The 
Vaisala method is useful for estimating an orbit from a short-arc set of observations spanning 
a day or two (Vaisala 1939). The technique assumes that the object is at apsis at either 
the starting or ending observation, so that the object’s instantaneous radial velocity is zero. 
It can provide a useful discriminator between Main Belt and NEA orbits when very few 
observations are in hand. Next, an attempt is made to hnd additional tracklets for that 
object by comparing residual positional differences between each orbit and tracklet pair 
(noting magnitude and consistency of astrometric residuals). Then, orbits are computed 
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using Gauss’ method (or other methods such as those described in, e.g., Marsden 1985) for 
linked tracklets. The process is repeated using a different starting orbit technique (c.f. the 
assumed elements or generalized Vaisala techniques discussed in Marsden 1991). At each 
step, confirmed linkages result in tracklets being removed from the tracklet stack, further 
simplifying linkage attempts going forward. 


Approximately 10,000 tracklets from the synthetic all-sky survey were evaluated in this 
fashion. In the procedure described above, all objects were treated as new and unidenti- 
hed. In practice, this is an unrealistic assumption since the current MBA catalog contains 
~500,000 objects with well-determined orbits. By the time a new space-based survey mission 
could be launched, even assuming construction started today, we expect the MBA catalog 
to contain ~1 million objects with well-determined orbits, and the NEO catalog could be 
expected to contain perhaps 15,000 - 20,000 objects. As a standard practice, the MFC runs 
all tracklets through procedures to check against known objects to remove in real time the 
easiest linkages to make. It is expected that >30% of tracklets will be immediately identifi¬ 
able, greatly simplifying the linking process. This result also serves as a warning that any 
dramatic increase in false tracklet generation rates will greatly complicate linking and thus 
reliability. It is imperative that the submitted tracklets be dominated by real objects. 


5. Conclusions 

We have constructed a set of simulations that account for updated NEA population 
models, including physical properties, orbital elements, and numbers. A portion of our 
model demonstrates the moving object pipeline’s ability to extract and link tracklets in the 
presence of potentially confusing astrophysical sources as well as other asteroids that are 
not NEAs. We have successfully computed accurate orbits for the assumed trial survey 
cadence using the standard tools employed by the MFC. However, we note that the ability 
to compute orbits depends critically on excluding most false detections from the submitted 
tracklets. Construction of a more complete catalog of potentially confusing non-NEAs such 
as MBAs and Jovian Trojans is essential. 

We compared the performance of the Venus-trailing and LI surveys and found that 
the Venus-trailing survey detects slightly fewer FHAs >140 m than the LI survey, even 
assuming that no degradation results from the inability to downlink full-frame images. While 
the Venus-trailing survey discovers more Amors than the LI survey, these objects are less 
likely to constitute impact hazards compared to Atens and Apollos, nor are they likely to be 
suitable targets for future missions. The LI survey discovers more lEOs, Atens, and Apollos 
than the Venus-trailing survey. These results demonstrate that the cost, complexity, and risk 








-42- 

associated with sending a survey telescope to a Venus-trailing orbit is unwarranted. While 
neither survey is capable of fulhlling the 2005 Congressional mandate to NASA to hud 90% 
of all near-Earth objects larger than 140 m in diameter by 2020, an advanced space-based 
survey can make signihcant progress quickly. 

No asteroid survey has demonstrated the ability to detect NEOs at low SNR and high 
reliability using only the equivalent of a spacecraft processor and memory to perform all 
source extraction algorithms. The loss of ancillary science and the ability to perform optimal 
image reprocessing with calibration products derived later in the mission associated with 
downlinking only regions of interest instead of full-frame images represents an additional 
penalty for interior-orbiting missions. 

These simulations serve as a starting point for optimizing the performance of an ad¬ 
vanced survey aiming to largely characterize the hazard posed by NEOs, identify the best 
targets for future exploration, and deliver a high-quality dataset ripe for solar system science 
to the community. Further improvements, such as those described above, can enhance the 
fidelity of the results. 
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